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Abstract 

This article surveys the development of observational understanding 
of the interior rotation of the Sun and its temporal variation over ap- 
proximately forty years, starting with the 1960s attempts to determine 
the solar core rotation from oblateness and proceeding through the devel- 
opment of helioseismology to the detailed modern picture of the internal 
rotation deduced from continuous helioseismic observations during solar 
cycle 23. After introducing some basic helioseismic concepts, it covers, 
in turn, the rotation of the core and radiative interior, the "tachocline" 
shear layer at the base of the convection zone, the differential rotation 
in the convection zone, the near-surface shear, the pattern of migrating 
zonal flows known as the torsional oscillation, and the possible temporal 
variations at the bottom of the convection zone. For each area, the article 
also briefly explores the relationship between observations and models. 
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1 Introduction 



The internal rotati on of the Sun is inti mately related to the processes that drive 
the activity cycle. iBrown et al. stated that, "Knowledge of the internal 



rotation of the Sun with latitude, radius, and time is essential for a complete 
understanding of the evolution and the present properties of the Sun," and this 
remains true today. 

The Sun rotates on its axis approximately once every twenty-seven days; 
however, the rotation is not uniform, being substantially slower near the poles 
than at the equator. This superficial aspect of the solar differential rotation 
was well known from sunspot observations as early as the seventeenth century. 
However it is only within the last thirty years that it has become possible to 
observe the rotation profile in the solar interior, and mostly within the most 
recent solar cycle that its subtle temporal variations have become evident. He- 
lioseismology - the study of the waves that propagate within the Sun and the 
inference from their properties of the solar interior structure and dynamics - is 
the most important tool we have to measure this internal rotation. 

In this review, we start by introducing some of the basic concepts of helioseis- 
mology (§[2]) and the inversion problem (§[3]) as it applies to the internal solar 
rotation. Next, after a brief historical overview (§ [4]) of the observations, we 
consider what we have learned from hclioscismology about the rotation profile 
and its variation with depth. 

We consider first the time-invariant part of the solar rotation profile. The 
main features of interest are (Figure [T]): 

1. the radiative interior and core, which appear to rotate approximately as 
a solid body, though the innermost core may behave differently; (§[5]), 

2. the tachocline, a relatively thin zone of shear between the differentially 
rotating convection zone and the radiative interior, which is believed to 
play an important role in the solar dynamo (§[6]); 

3. the differential rotation in the bulk of the convection zone (§[7]); and 

4. the subsurface shear layer between the fastest-rotating layer at about 
O.95i?0 and the surface. (§[8l). 

We will consider each of these in turn, working outwards from the core to the 
surface, and then discuss the time-varying part of the rotation - the torsional 
oscillation (§ ^ and the possible variations at the base of the convection zone 
(^fTO|. We attempt to place the observations in the context of models; however, 
this is a review from an observer's point of view, and an exhaustive examination 
of the models themselves is beyond its scope. 
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Figure 1: A section through the interior of the Sun, showing the contours of 
constant rotation and the major features of the rotation profile, for a temporal 
average over about twelve years of MDI data. The cross-hatched areas indicate 
the regions in which it is difficult or impossible to obtain reliable inversion results 
with the available data. 
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Figure 2: A single Doppler velocity image of the Sun from one GONG [Global 
Oscillation Network Group] instrument (left), and the difference between that 
image and one taken a minute earlier (right) . with red corresponding to motion 
away from, and blue to motion towards, the observer. The shading across the 
first image comes from the solar rotation. 



2 Acoustic Modes 

The raw data of helioseismology consist of measurements of the photospheric 
Doppler velocity - or in some cases intensity in a particular wavelength band 
- taken at a cadence of about one minute and generally collected with as little 
interruption as possible over periods of months or years; the measurements can 
be either imaged or integrat ed ( "Sun a s a Star" ) . An overview of the observation 
techniques can be found in iHill et al\ ([l99l3). Figure [5] shows a typical single 



Doppler velocity image of the Sun, and Figure[3]a portion of an Z = time series, 
derived by averaging the velocity over the visible disk for each successive image 
in a set of observations. The five-minute period and the rich beat structure 
are clearly visible in the time series. For an example of an integrated-sunlight 
spectrum from a long series of observations, see Figure fTSl 

As was first discovered bv lDeubner (1973), the velocity or intensity variations 



at the solar surface have a spectrum m k — lu or I — v space that reveals their 
origin in acoustic modes propagating in a cavity bounded above by the solar 
surface and below by the wavelength-dependent depth at which the waves are 
refracted back towards the surface. These "p modes" can be classified by their 
radial order n, spheric al harmonic degree / , an d azimuthal order m; as discussed, 
for example, in § 2.2 of Aaron C. Birch ( 20051) . the radial displacement of a fluid 
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Figure 3: A segment of an ^ = time series of Doppler velocity observations, 
showing the dominant five-minute period and the rich beat structure. 
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element at time t, latitude 9 and longitude (f> can be written in the form 

I 

6r{r,9,cf>,t)= ^ a,uUnl{r)Yr{eA)e'^""-\ (1) 

m— — l 

where ^nim is the radial eigenfunction of the mode with frequency tOnim and 
Yi"^6, (/) is a spherical harmonic. As seen in Figure HI the power in the spectrum 
falls along distinct "ridges" in the I — v plane, each ridge corresponding to 
one radial order. The modes making up the n = ridge are the so-called / 
modes, which are surface gravity waves. The p modes, so called because their 
restoring force is pressure, are excited at the surface and have their largest 
amplitudes there. Another class of modes, the g modes with gravity as the 
restoring force, excited in the core and with amplitudes vanishing at the surface, 
are hypothesized to exist but have so far not been definitely observed (ij 15. 9p . 




Figure 4: Typical I — v spectrum from one day of GONG observations. (Image 
courtesy NSO/GONG.) 

The longer the horizontal wavelength - and the lower the degree - the more 
deeply the modes penetrate, with the radial I — mode going all the way 
to the core of the Sun (but providing no rotational information), while modes 
with I > 200 or so penetrate only a few megameters below the surface and 
are generally too short-lived to form global standing waves; these are the modes 
used for local helioseismology. The lower turning point radius, rj, is a monotonic 
function of the phase speed v/L, where L — + 1) ~ I + 1/2, as shown in 
FigureO The varying penetration depth with degree, as illustrated in Figure[6l 
makes it possible to deduce the rotation and other properties of the solar interior 
profile as a function of depth. 

The lowest-degree modes are observed in integrated sunlight, but for the 
purposes of measuring the interior rotation profile we are mostly concerned 
with what are termed medium-degree (l < 300) modes, which can be observed 
with imaging instruments of relatively modest (« 10 arcsec) resolution. The 
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Figure 5: Lower turning point of acoustic modes as a f unctio n of phase speed 
vjL, based on Model S of Christensen-Dalsgaard et al. ( 1996f ). 



power in the modes peaks at about 3 mHz, or a period of 5 minutes; useful 
measurements can be made for modes between about 1.5 and 5 mHz, with 
the frequency determination becoming more challenging at the extremes due to 
signal-to-noise issues and, at the high-frequency end, to the increasing breadth 
of the peaks. 

2.1 Differential rotation and rotational splitting 

The Sun's rotation lifts the degeneracy between modes of the same I and different 
TO, resulting in "rotational splitting" of the frequencies as waves propagating 
with and against the direction of rotation (prograde and retrograde) have higher 
and lower frequencies. To first order, the splitting 5i'rn,i = i^-m.i ~ t^+m.i is 
proportional to the rotation rate multiplied by to. 

Figure [7] shows a typical m—v acoustic spectrum of GONG data at I — 100. 
The effect of the rotation causes the ridges at each n to slant away from the 
u = Q axis; closer examination reveals that the ridges have an S-curve shape 
that arises from the differential rotation, and also shows the ridge structure due 
to leakage, which will be discussed below in § 12.21 

Because modes of different to values sample different latitude ranges, with 
the sectoral (|to| — I) modes confined to a belt around the equator and the 
zonal or TO = modes reaching to the poles, as illustrated in Figure [8l we can 
measure the rotation as a function of latitude. 

A given {n, I) multiplet consists of 2Z -f 1 frequency measurements if each 
{I, to) spectrum is analyzed separately, though some fraction of these frequencies 
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Figure 6: Locations of modes in the I, v plane for a typical MDI mode set. The 
colored shading represents the radial regions in which the modes have their lower 
turning points; the core, r < Q.2Rq, the radiative interior, 0.2 < r/RQ < 0.65, 
the tachocline, 0.65 < t/Rq < 0.75, the bulk of the convection zone, 0.75 < 
t/Rq < 0.95, and the near-surface shear layer, t/Rq > 0.95; the dashed line on 
the lower right corresponds to r/RQ = 0.99. 
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Figure 7: Spectrum for I — 100 in the ^, m plane (top) and detail (bottom) of 
a single ridge (radial order) showing the curvature due to differential rotation 
and the multiplc-ridgc structure arising from spherical harmonic leakage. 




Figure 8: Spherical harmonic patterns for I — 50: (left, to = 0; center, m — 45, 
right; m = 50). 
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may be missing in any given data set. This amount of data was somewhat 
unwieldy in the early days of helioseismology. It is therefore common to express 
I'niinT') as a polynomial expansion, for example, 



l^nlm = ^nl + > % (n, /) P ■ ^m) , (2) 



where the basis functions are polynomials related to the Clebsch-Gordan coef- 



ficients C'^^ by 



(jRitzwoller and Lavelvl . Il991r ). Indeed, in many analysis schemes coefficients of 
the expansion are derived by fitting directly to the acoustic spectrum and the in- 
dividual frequencies are not measured. This approach can improve the stability 
of the fits, perhaps at the cost of imposing systematic errors. Early work used 
Legendre polynomials; however, most modern work uses either Clebsch-Gordan 
coefficients or the RitzwoUer-Lavely formulation, which come closer to being 
truly orthogonal for the solar rotation problem. Only the odd-order coefficients 
encode the rotational asymmetry, while the even-order coefficients contain in- 
formation about the structural asphericity. Roughly speaking, the ai coefficient 
describes the rotation rate averaged over all latitudes, and the 03 and higher 
coefficients describe the differential rotation. 



2.2 Spherical Harmonics and Leakage 

Spherical harmonic masks are used to separate the contributions from modes of 
different degree and azimuthal order into complex time series, which can then 
be transformed to acoustic Fourier spectra. 

The radial component of the velocity at the solar surface from a mode with 
a given degree /, azimuthal order m and radial order n is given by 

Vn,i,m{<jy,0,t) = Re[a,u,m(t)Pz''"'(cos0)e™^)], (4) 

whe re ReH denotes the real pa rt, is longitude and 6 is latitude. (See, for exam- 
ple, Schou and Browiil 1994al ) The masks used separate the different spherical 



harmonics take the form 



Mum'^Yi^ra{eA)A{p), (5) 

where A is an apodization function and p = \/ cos^ d + sin^ d sin^ </> represents 
the fractional distance from disk center in the solar image. The line-of-sight 
projection factor \&V — \/l — f?. 

Because only part of the solar surface is visible at any time, the masks are 
not completely orthogonal and the modes "leak" into neighboring spectra. This 
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leakage complicates the analysis and can cause systematic errors in the measured 
frequencies if it is not correctly taken into account. For a detailed discussion o f 
the calculation of the so-called "leakage matrix," see Schou and Browiil ( 1994a ): 



Hill and Howe ( 1998f) . Briefly, the leakage matrix element s(/, m, m') for leak- 



age from the Z', m! mode to the Z, m spectrum can be computed as 



1 



1 ^7r/2 



s{l,m,l',m') ^ - / Pl''{x)PP' {x)cos{m(t))cos{m'(j))V{p)A{p)dxd(j). 



I" J-l J-7r/2 

(6) 

Symmetry properties in this expression lead to some simple exclusion rules; for 
example, leaks with odd \Sl + Sm\ (where Sm = m ~ m' and SI = I — I') are not 
allowed. 

One example of the importance of the leakage is in the contribution of the 
so-called m-leaks {61 = 0,6m = ±2) to the estimation o f low- degree splittings. 
As pointed out, for example, bv iHowe and Thompson! ( 1998f ). these leaks are 



strongest for small \m\; they are also asymmetrical, especially for \m\ = I, where 
the m = I peak has an m = I — 2 leak on one side and no counterbalancing 
m = I + 2 leak on the other. Especially for 1 = 1, this can introduce a serious 
systematic error into the estimate of the splitting if not properly accounted for. 

Leakage also means that integrated-sunlight instruments (which effectively 
use the / = mask) can detect modes of < Z < 5, though the sensitivity falls 
off rapidly for / > 1. All these modes appear in a single acoustic spectrum; the 
instruments are sensitive to odd-m modes for odd I and to even-m modes for 
even I, with the sectoral, or \m\ = I, modes most strongly detected. 

In general, the leakage has effects throughout the acoustic spectrum, but the 
most deleterious effects arise when the leaks cannot be resolved from the target 
peaks. This occurs for m-leaks at frequencies above about 2 mHz; for higher- 
degree modes the leakage between modes of adjacent I becomes a problem, as 
the ridges become both broader, and more closely spaced in frequency, at around 
I = 150. Beyond this point the peaks cannot be fitted independently, and some 
modeling of the leakage is essential in order to estimate the mode parameters. 

2.3 Estimating rotation properties directly from coeffi- 
cients 

It is possible to make some inferences about the rotation profile without carrying 
out a full-scale inversion. Simple examination of the odd-order coefhcients, 
sorted by the lower turning-point radius of the modes, reveals the existence 
of the near-surface shear, the differential rotation within the convection zone, 
and a discontinuity in the differential rotation at the base of the convection 
zone, as shown in Figure [9l More soph isticated analysis is also possible. For 



example, IWilson and Burtonclav ( 1995 ) gave approximate expressions for the 



rotation profile at different latitudes as sensed by a particular n, I multiplet, 
n™', as follows: 

^af +af +af, (7) 
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Figure 9: ai (top) and 03 (bottom) coefficients for (left), 1986 BBSO observa- 
tions, (middle) 108 days of GONG observations in 1996, (right) the mean of 35 
consecutive 108-day periods of GONG observations from 1995-2005, plotted as 
a function of phase speed with the turning point radius marked on the upper 
axis. 
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These estimates, where the subscripts on the L HS refer to the latitude in de - 
grees, are noisy for individual multiplets, but Wilson and Burtonclav ( 19951 ) 
were able to build up a picture of the internal rotation from BBSO data by 
forming cumulative averages with the input data sorted in ascending order of 
v/L. 



13 



3 Inversion Basics 



Various inversion techniques are used to infer the internal rotation profile from 
the observed frequency splittings. The aim of the inversion procedure is to 
form linear combinations of the data that give well-localized inferences of the 
rotation at a particular location within the Sun. We will discuss only linear 
inversion methods, as non-linear approaches are not needed for the relatively 
low velocities involved in the global rotation 



3.1 The inversion problem 



The basic 2-dimensional rotation inversion problem can be stated as follows: we 
have a number M of observations d,;, from which we wish to infer the rotation 
profile f2(r, 9) where r is distance from the center of the Sun, and 9 is (conven- 
tionally) colatitude. Each datum is a spatially weighted average of the rotation 
rate: 



Ki{r, 9)n{r,9)drd9 



(11) 



where Rq is the solar radius, the error term e corresponds to the noise and 
measurement error in the dat a, and K is a model-dependent spa tial weighting 
function known as the kernel (jHansen et al\ . 11977 ; ICuvperd . lT980l) . For the two- 
dimensional rotation inversion, the radial part is related to the eigenfunction 
of the mode and the latitudinal part to the associated Legendre polynomial; 
Schou et ali (| 19941 ) give the expression for the kernel as 



Knlm{r,9) 



Prix)' 



(12) 



where 



dPT 

dx 



(1 - x') - 2P[ 



dx 



-X + 



1 — x- 



■prixf 



p{r)rsin9, 



Jo 



+ ?7„i(r)^]p(r)r^dr. 



(13) 



X — cos9, = l{l + 1), is the radial displacement for the eigenfunction of 
the mode, L^^rjni is the horizontal displacement, and p{r) is the density. (See 
Figure fTOl for illustrations of sample kernels.) 
The aim of the inversion is to find 



M 



n{ro, 9o) = ^ Ci{ro, 9o)di 



(14) 



where (rg, 9q) is the location at which the inferred rotation rate fl is to be found 
and the Ci are the coefRcients to be used to weight the data; the inversion process 
can be thought of as the search for the best values for these coefficients. 
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Figure 10: Sections through rotation kernels for selected azimuthal orders for 
I = 3,n = 9 (top) and I ~ 20, n = 5 (bottom). 
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3.2 Averaging kernels 

By substituting Equation pT|) into the RHS of Equation (THl) we obtain 

^^(^0,^0)= / ° I IC{ro,9o;r,e)n{r,9)drd0 + e^, (15) 







where 



M 



IC{ro,9o; r, 0) = J2 c^iro,9o)K,{r, 6) (16) 



is the averaging kernel for the location (ro,0o)- The averaging kernels are in- 
dependent of the values of the data. However, because the uncertainties in the 
data are used to weight the inversion calculation that generates the coefficients 
Ci, as described below in 13.41 and l3.5[ these do enter indirectly into the aver- 
aging kernels. The averaging kernels also depend on which modes are present in 
the input data set. They provide a useful tool for assessing the rel iability of an 



inversion in ference from a particular mode set (see, for example, ISchou et al. 



19921119941) 



3.3 Inversion errors 

If the errors on the input data are uncorrelated and properly described by a 
normal distribution whose width corresponds to the quoted uncertainty cr^, the 
formal uncertainty on the inferred profile is given by 

a^Mro, 9o)] = ^[c,(ro, 0o)^t^]'. (17) 

i 

In the (usually unrealistic) case where the errors on the input data are all equal, 
we can write 

a^[niro,9o)]^Airo,9o)a, (18) 
where the "error magnification" is given by 

A{ro,9o)^J2^c.,iro,9of]'^'. (19) 



As discussed, for example, bv lChristensen-Dalsgaard et al. ( 1990h . a quantita- 



tive choice of regularization parameters can then be made by finding the "knee" 
of a tradeoff curve where the error magnification is plotted against the width of 
the averaging kernel. However, in the two-dimensional case this does not always 
give a clear result, and this formulation of the error magnification is not very 
useful for modern data sets where the the uncertainties on the parameters are 
anything but uniform. Instead, one can consider the uncertainty on the inferred 
quantity at a particular location. 

Even when the errors on the input data are uncorrelated, the errors on the 
inferred profile will not be, as discussed by Howe and Thompson ( 19961 ). (As 
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a simple way to understand this, consider the case where one measurement is 
significantly "off"; this will affect the inferred profile at every location where 
the inversion coefficient Ci for that datum is non-zero.) In the one-dimensional 
case, the correlation between the errors for two points tq and ri is given by 

r(r. ) = Ec»(ro)c,(ri)af 



this ca n easily be generalized to the two-dimensional case. iHowe and Thompson 



(Il996h Found that the spatial scale over which the inversion errors are signifi 



cantly correlated is usually similar to that for the averaging kernels, though for 
some cases where the inversion parameters have been badly chosen the results 
can be correlated over long distances even when the averaging kernels appear 
well formed. 

Error correlations by definition should not distort the inferred profile beyond 
the distribution predicted by the formal uncertainty on the inferences, provided 
always that the input uncertainties are correct. However, the finite width of 
averaging kernels also gives rise to a systematic error that can be much larger. 
Consider, for example, the case where a thin shear layer is not resolved; then all 
the estimated rotation rates on one side of the shear could be underestimated, 
and those on the other side overestimated, by several times the formal uncer- 
tainty. Such systematic errors and t heir relationship to the averaging k ernels 
ha ve been discussed, fo r example, by Christensen-Dalsgaard et al. ( 199(J ). 



Gough et al. ( 1996[ ) pointed out that it is not sufficient for the rotation rates 



at two locations to have non-overlapping errors as calculated in Equation (|17p . 
and described a method for increasing the error estimates on inversions to al- 
low truly significant differences between the inferred rotation rate at different 
locations to be determined. This method, however, has not been widely used. 

Because the input data are noisy and of finite resolution, the inversion prob- 
lem does not have a unique solution; there will always be a tradeoff between 
noise and good localization. Two widely-used approaches to balancing these cri- 
teria are "regularized least squares" (RLS) and "optimally localized averaging" 
(OLA). 

3.4 Regularized least squares 

The RLS approach to the inversion problem is to find (essentially through a 
least-squares fit) the model profile that best fits the data, subject to a smooth- 
ness penalty term, or regularization. More regularization - a larger weight- 
ing for the penalty term - results in poorer spatial resolution (and potentially 
more systematic er ror) but smaller uncertainties. In one such implementation 
( Schou et al\ . Il994l ) , we minimize 



(21) 
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with /ir and ne being the radial and latitudinal tradeoff parameters. The RLS 
inversion has the advantages of being computationally inexpensive and always 
(thanks to the second-derivative regularization, which amounts to an a priori 
assumption of smoothness) providing some kind of estimate of the quantity of 
interest even in locations that are not, strictly speaking, resolved by the data. 
In this method, the averaging kernels /C can (but need not be) calculated from 
the coefficients in a separate step. They are not guaranteed to be well localized, 
though they are forced to have a center of mass at the specified location rg, ^q. 
Figure [TT] illustrates typical averaging kernels for a 2dRLS inversion of an MDI 
data set. 



3.5 Optimally localized averaging 

In the Subtractive OLA (SOLA) approach (jBackus and Gilbertl . [l96lll970( ). the 
minimization is applied to the difference between the actual averaging kernels /C 
and a target kernel T, for example a 2-dime nsional Gaussian or Lorentzian func- 
tion. In this case ( Piipers and Thompsonl . 1992 . .1994) the function minimized 



^rdrdO + [o'iCi{rQ, 
1=1 



d?- (22) 



Both the tradeoff parameter A and the radial and latitudinal resolution of the 
inversions must be chosen before running the inversion. If the choice of target 
kernel is poor - too narrow or too wide for the quantity and quality of the data 
- the reliability of the inversion will suffer. In OLA inversions, setting target 
locations outside the regions that can be resolved using the data will result in 
averaging kernels displaced from their targets, and this should be taken into 
account when interpreting the results. Figure [T^ illustrates typical averaging 
kernels for a 2dS0LA inversion of an MDI data set. 

Another approach, older, and m ore computationally expensive, is the Multi- 
phcative OLA (MOLA) described bv lPiipers and Thompsonl (|l992lll994l ). Here, 
no target form is imposed on the averaging kernel, but it is multiplied by a term 
which penalizes large values away from the target location. 



3.6 Other methods 

Alternatives to full 2-dimensional inversions are the so-called "1.5-dimensional" 
approach, in which 1-dimensional radial inversions are carried out separately 
for each of the coefficients describing the latitudinal rotation variation, and 
"IdiE" Id" inversions in which the radial and latitudinal variations of the rotation 
rate are integrated separately. For details of many of these methods, please see 
Schou et and references therein. 
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Figure 11: Averaging kernels for a typical RLS inversion of MDI data, for 
target latitudes (a), 15 (b), 30 (c), 45]^^), 60 (e) and 75 (f) degrees as marked 
by the dashed radial hnes, and target radii 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 
0.99 Rq indicated by colors from blue to red as denoted by the dashed concentric 
circles. Contour intervals are 5% of the local maximum value closest to the 
target location, with dashed contours indicating negative values. 




Figure 12: As Figure [TT| for a SOLA inversion. 
20 



3.7 Limitations 



It is important to bear in mind the limitations of the inversion process when 
considering the results. The deepest and shallowest depths that can be resolved, 
for example, are limited by the deepest and shallowest turning-point radii of 
the available modes. The rotational splitting at a given m is to first order 
proportional to the rotation rate multiplied by \m\; since the only mode whose 
latitudinal kernel reaches the pole is the m = mode, which has no longitudinal 
structure and so can convey no rotational information, and the modes of small 
\m\/l have only a few nodes around the equator and hence have low sensitivity 
to the rotation, the 2d inversion becomes progressively less reliable at high 
latitudes. Furthermore, since only modes of relatively low degree {I < 20) 
penetrate into the radiative interior, the latitudinal resolution in this region is 
quite poor and becomes progressively worse with depth; radial resolution also 
becomes coarser in the interior. The practical effects of such limitations can be 
assessed by careful inspection of the averaging kernels, or by performing forward- 
calculation tests in which the averaging kernels are convolved with known test 
profiles. 

Another point to bear in mind when considering inversion results is that the 
inversion can measure only the north - south symmetric part of the profile; any 
asymmetry between the hemispheres is averaged out. The inversions are also 
insensitive to meridional motions. Some information on hemispheric differences 
cal l be obtained us i ng th e the techniques of local helioseismology, as reviewed 
by lAaron C. BirchI ( 2005h . but these techniques, using high-degree modes, are 
mostly sensitive only to the outer layers of the Sun. 
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4 Observations: A Brief Historical Overview 

Systematic helioseismic observations stretch back nearly thirty years, as illus- 
trated in the schematic chart in Figure 1131 
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Figure 13: Schematic time line of helioseismic observations in the last three 
solar cycles (top panel), with the filled part of each bar representing approximate 
duty cycle, plotted on the same temporal scale as the butterfly diagram (bottom 
panel) of the gross magnetic field strength from Kitt Peak observations. 



Prior to the identification of global low-degree modes by IClaverie et 



( 1979f) . observing runs were usually short and carried out at a single site 
However, the advantages of more extended observations (to obtain better fre- 
quency resolution), and of observations no t mo dulated by the day-night cy - 
cle, were soon r e cogni z ed. Grec et al. ( 1980[ l and Duvall Jr and Harvev (198J); 
Duvall Jr et al\ ( 1984L 1986[ l carried out important observations at the South 
Pole during the Austral summer, but for long time series it is more practical to 
observe either from a network of sites spaced around the world, or from space. 

Some of the first long-term sets of low-degree observ ations came from the 
Active Cav ity Irradiance Monitor (ACRIM) experiment ( Woodard and NoyesL 
aboard the Solar Maximum Mission spacecraft, which took helioseis- 
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mic measurements in 1980 and 1984-1985, the Mark I instrument in Tenerife 
( Palle et al ]. ll989l). and the pre c ursors of the Birmingham Integrated Solar Net- 



work (BiSON) (|Elsworth a/.1 . Il990al) . Meanwhile, resolved-Sun observations 



were carried out at the South Pole by Duvall and collaborators, and by various 
other observers in the USA; these observations will be discussed in more detail 

latere 

Libbrecht and WoodardI (|l990[ ) observed the medium-degree modes from 



Big Bear Solar Observatory (BBSO) in the 1986-1990 rising phase of so- 
lar cycle 22. The first observations from widel y separated sit e s wer e carried 



out by the Birmingham/Tenerife group in 1981 (|Claverie et al\ . \l9S4 ). and by 



1992 the six-station BiSON network was complete; it has been operating ever 
since. Another ne twork of integrated-sunlight instruments, the French-based 
IRIS (|Fossati[l995l) . operated from 1989-2003. 

The Global Oscillation Network Group (GONG) (|Harvev et all Il99d ) has 
been collecting continuous, high-duty-cycle observations of the medium-degree 
p modes since 1995, using a six-stati on worldwide network , and the Michelson 



Doppler Imager (MDI) instrument (jScherrer et all Il995f ) aboard the SOHO 



spacecraft has been in operation since 1996, so that these two projects have 
essentially complete coverage of solar cycle 23. SOHO also carries instruments 
dedicate d to the study of low -degree oscillations; LOI (Luminosity Oscillations 
Imag er) (Frohlich et al . l995l ). and GOLF (Global Oscillations at Low Frequen- 
cies) ( Gabriel et all Il995l ). This wealth of high-quality data has given us the 
opportunity to study the solar interior rotation and its solar-cycle changes in 
more detail than ever before. 

T99I 



Also worth noting are the LOWL-ECHO project (jTomczvk et 



which made medium-degree observations from one or two sites from 1994 to 
2004. and the high-degree Taiwanese Oscillations Network ( Chou et al. . 19951) 
deployed over the 1993-1996 period. 

All these observations will be considered in more detail as we proceed to 
examine the results pertaining to the interior rotation. 
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5 The Core and Radiative Interior 



5.1 The oblateness controversy 



Interest in the rotation of the deep solar interior predates systematic hehoseismic 
observation. One other possible diagnostic of the internal rotation is provided 
by the solar oblateness; because the Sun is not a solid body, both gravitational 
and rotational effects cause a very slight flattening. The lowest-order term in 
this effect is related to the quadrupole moment J2; confusingly, the next-highest 
te rm, J4, is sometimes called octopole and sometimes hexadecapole. According 
to iRozelot and RoeschI (|l997r ) , who give a useful review of attempts to measure 
the solar oblateness, for a non-rotating Sun the oblateness Ar = r^g — Vp^i is 
given by 

- - (23) 

where Veq and rpoi are the equatorial and polar radii, respectively, and ro is the 
radius of the best sphere passing through req and Vpoi . If there is an additional 
6r contribution from the surface rotation this expression becomes 



Ar — Sr 



I*. 

ro 2 



(24) 



Th e units of 5r and Ar are conventionally arc ms. 

Dickd (1964) noted that, if the Sun were oblate because of fast interior ro- 



tation, the effect on its gravitational potential might destroy the agreement 
between the predictions of General Relativity and the observations of the per- 
ihelia of the inner planets, (specifically Mercury, though Venus could in princi- 
ple experience a smaller effect) potentially leaving room for alternate theories 
of gravitation. Dicke set out to determine the solar oblateness from ground- 
based measurement s - a challeng ing endeavor that produced controversial re- 
sults. Models fe.g.. iBrandtl 119661 ) suggested that the interior of the Sun could 
still be spinning at the rapid rate at which it originally formed, while the exte- 
rior had been slowed down by the torque of the solar wind. (As will be further 
discussed in Section [51 in the absence of direct observations of the solar interior 
the picture of solar interior dynamics was not at all clear, although the exis- 
tence of something like what we now call the tachoclin e could be inferred from 
theoretical arguments.) lOicke and Goldenberd ( 1967b[ ) reported finding a solar 
oblateness value of 5 x 10~^, which would be sufficient to create an 8% discrep- 
ancy between observations and the Einsteinian prediction for the precession of 
the perihelion of Mercury, and would imply a fast-rotating core. 

The results, and the inferences Dicke and collaborators drew from them, 
raised a storm of controversy that may well have helped to stimulate interest in 
the Sun's interior rotation profile. The criticisms and Dicke's responses to them 
would fill a lengthy rev iew article by themselves; we give only a few examples 
here. iRoxburghl (|l967l ) suggested that the re sult might be explaine d by th e 
solar different i al rot ation, an idea rejected by iDicke and Goldenberg (1967a). 
Howard et ali (|l967n concluded, on the basis of a variety of simple models of the 
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solar "spin-down," that the Sun should h ave reached a state of uniform rotation 
quite quickly after its initial formation. ISturrock and Gilvarrvl (jl967l ) pointed 
out that the presence of magnetic field in the s olar interior might wel l com - 
plicate the issue, and in an accompanying article iGilvarrv and Sturrockl (|l967l ) 
suggested using a space probe in a highly eccentric orbit as a more direct test 
of general relativity - or, alternatively, that "more complete theoretical and 
observational knowledge of the visible layers and the interior of the Sun" was 
needed. 

At least partly inspired by the controversy, iKraftI (jl967l ) studied the rota- 
tional velocities of young solar-type stars in the Pleiades and concluded that 
angular momentum was lost on a timescale of about half a billion years, but 
noted in his conclusion that "it is w rong to conclude that the pres ent work in 
any way supports the Dicke result." iGoldreich and SchubertI (|l968l ) considered 
the stability of differentially rotating stars and concluded that it was possi- 
ble but not likely that a radia l rotation gradient such as that required by the 
Dicke and Goldenbers (|l967bh result might exist. 



H. 



Hill, a former colleague of Dicke who had he l ped b uild the instrument 



nad ne ip 

with which the 1964 observations were made ( Dickd . fl964l ). and collaborators, 
also attempted to measure the solar oblateness, using an instrument, SCLERA 
[Santa Catalina Observatory for Experimental Relativity by Astrometry] , which 
was later to play a r ole in the early days of he lioseismology. This measurement, 



carried out in 1973, ()Hill and Stebbinsl■ll975^■ found a 9.6 X 10 ^ value for the 



oblateness, much smaller than that of lDicke and Goldenberg|(|l967bh : [Hiire^ al 



( 19741 ) also pointed out a time- varying difference between the brightness of the 
solar limb and poles that might account for the anomalously high oblateness 

me asurement. 

Ulrich and Hawkind ()l981al lbh made an early attempt to deduce what the J2 



and J4 terms should be based on a simple differential rotation profile deduced 
from surface measurements, obtaining predicted values of between 1 and 1.5 x 
10"'' for J2 and between 2 and 5 x 10^^ for J4 depending on the size of the 

co nvective enve l ope. 

Dicke et 'd\ (Il986l . ll987l) repeated the 1966 measurements with an improved 



instrument, and obtained significantly smaller values for the obl a tenes s, with 
some weak evidence for a solar-cycle variation. Lvdon and Sofia ( 1996f ) made 
measurements using a balloon-based instrument and obtained values of 1.8 x 
10"'' for J2 and 9.8 x 10^'' for J4. By this point, however, the focus in the 
solar oblatene s s stud ies had moved away from trying to infer the core rotation. 



Mecheri et al\ (j2004[ ) used more realistic models of the internal rotation profile to 



suggest that the J4 term should be particularly sensitive to the subsurface shear. 
Recent work on determining the oblateness from the shape of the solar limb 
has taken i nto account considerations of ne a r-surf ace temperature or magnetic 
variations. iKuhn et al\ (|l998r ). lEmilio et al\ (120071) used obse rvations from MDI 
during rolls of the SOHO spacecraft and lFivian et al\ (|2008l ) used the RHESSI 
X-ray telescope. The work with SOHO revealed a temporal variation in the 
shape of the solar limb, with greater apparent oblateness at solar maximum, 
suggesting that hotter, brighter activity belts have greater apparent diameter. 
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This poses an apparent contradiction to the results obtai ned from hehoseismic 
inferences of the asphericity. Indeed, iFivian et all (|2008f ) suggest that all the 
temporally- varying, excess oblateness found in the observations can be corrected 
away by removing an ad-hoc term related to magnetic elements in the enhanced 
network. 

Meanwhile, a much more flexible tool - helioseismology - had become avail- 
able for probing the interior solar rotation. 



5.2 Early low-degree helioseismic results 

Around the early 1970s there were numerous attempts to search for global p- 
mode oscillations, with interest at first focusing on longer-period oscillations, the 
low-order, low-degree modes. Various theoretical predi c tions (jScufiaire et al 



llben Jij . Il976t IChristensen-Dalsgaard and Goughl . Il976l) of the periods 
were available, offering the hope that global oscillations could be used to probe 
the rotation and structure deep inside the Sun . At first most of the results 



( Livingston et al . 1977; Musman and Nve . 1977t Grec and Fossat 



negati ve, except for the 1 60- minute period of Severnvi et a/.l ( 1976f l: lBrookes et - 
( 19761 ). which was later ( Elsworth et al. . 19891 ) determined to be spurious and 



19771) 



were 



will not be further discussed here. The SCL ERA group ([Brown et all . 11978 



Hill and Caudell . I1979I : ICaudell and Hill . 1980l ) found a variety of longer-period 
fluctuations in their s olar-diameter data, b ut these results were not universally 
accepted; for example iFossat et ali (|l981al see also references therein) claimed 
that the SCLERA results were consistent with pure noise. 

Low-degree helios eismology became a reality when the Birmingham group 
(jClaverie et a/.l . ll979l ) identified oscillations in the five-minute frequency band in 
integrated sunlight as low-degree global modes, using observations from Tenerife 
and Pic du Midi during the summers of 1976-1978; these initial data were 
adequate only to identify the spacing between modes of the same I and different 
n, without resolving separate I = and I = 1 peaks. 



A French- American team ( Grec et aL . 11980*: 'Fos sat et al. . 1981b ) obtained 



five days of continuous observations at the South Pole in the austral summer 
of 1979-1980, and were able to identify peaks of degree 0, 1, 2, and 3 and 
even a weak / = 4 peak by superposing sections of the acoustic spectrum with 
different radial order. These modes were identified as being of radial order 
around 12-30, as opposed to the very low-order modes that had been sought 
in the low-frequency spectrum; both the noise characteristics of the spectrum 
and the low amplitude of the lower-order modes mean that the fundamental 
(/ = 0, n = 0) mode remains unobserved to thi s day, although some l ow- degree 



1996bl ). 



modes with single-digit n have been identified (Chaplin et al. 

Soon, the Birmingham team ( Claverie et ali 1981) . using 28 days of integrated- 
sunlight data from the Tenerife site and an analysis that involved "collapsing" 
segments of the acoustic spectrum so as to average together modes of the same 
degree and different radial order, reported finding three rotationally split com- 
ponents in the I = 1 modes and five in I = 2, with an average separation of 
0.75 /iHz. If correct, this would have implied a solar core rotation substantially 
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faster than the surface. Ilsaaki (|l982l ) suggested that the excess component peaks 
(when two and three would be the expected number for I — 1 and I = 2 respec- 
tively) could be explained if the solar core were rotating on an oblique axis^and 
had a very strong magnetic field; this idea, which was also mooted by iDicl^ 
( 1983f ) to explain a n oscillation of about half the solar rotation period seen i n 
the oblateness data (iDicke . 1979), was rebutted in some detail bv lGoughl(|l982l) . 

|l981bl ) reported that initial results from 5 days of low-degree 



Fossat et al 



observations at the South Pole suggested quite short lifetimes, about 2 days; the 
/ = peaks appeared narrower than those of / = 1 and I = 2. iGrec et al\ (Il983[ ) 
later identif ied about 80 normal modes in the South Pole data, but did not 
confirm the IClaverie et all (|l98ll ) rotational splitting result, instead reporting 
that the I = 1 peak seem ed too narrow to accommodate the reported splitting. 



Claverie et ali (|1982I ) reported a periodicity of approximately 13 days in the 



radial solar velocity, as measured using the resonant-scattering technique and 
the potassium D-line, and interpreted this as an e ffect of the solar core ro t ation ; 
however, this effect wa.s quickly explained away ("Purrant and Schroetei . 1985 



Andersen and Maltbvl . ll983t Edmunds and Gough . J983i : iDuvall Jr et al. 



as an artifact caused by the rotation of surface features - sunspots and plage - 
across the disk. 

Meanwhile, the low-degree five-minute acoustic spectrum had also been ob- 
served using the Active Cavity Irradi ance Monitor (ACRIM) aboa r d the Solar 



Maximum Mission (S MM) spacecr a ft (IWoodard and Hudsonl . ll983al ). IWoodard and Hudson 
( 1983bl ) agreed with IFossat et all ( 1981b ) in finding that the modes had hfe- 
times of abou t two days, too short for the rotational splitting reported by 



Claverie et all (Il98ll) to b e real. 

Later work (|Libbrechtl . Il988al : lElsworth et~al\ . Il990bl : IChaphn et aH [iM^) 
revealed that the width of the peaks - inversely proportional to the mode life- 
times - was strongly dependent on frequency across the five-minute spectrum, 
with lifetimes of a few days in the middle of the five-minute band and weeks or 
months at low frequencies where, unfortunately, the amplitudes of the modes 
are also small. Reliable direct measurement of the low-degree splittings would 
have to wait for some years, while sufficiently long, high-quality time series of 
data accumulated. 



5.3 Resolved-Sun measurements 



In the meantime, resolved-Sun observa tions provided some in f ormat ion about 
the rotation in the radiative interior. iDuvall Jr and HarvevI (|l984l ) reported 
observations at Kitt Peak, from 10-26 May 1983, for degrees 0<l< 100. When 
plotted as a function of degree, the results show a slow decrease in the rotational 
splitting, from the highest degrees down to about / — 6, with an unexplained 
bump at ^ = 1 1 , followed by an incre ase at lower degr e es up to a value of 660 nHz 
for Z = 1. These data, inverted bv iDuvall Jr et cd. (1984), yielded a rotation 
profile with much of the radiative interior rotating at or below the surface rate, 
but with a modest increase in the interior. A similar pattern was found by 



Browij (jl985f ). using 6 days of observations from the newly-developed Fourier 
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Tachometer, a true 2-dimensional imaging instrument that gave access to ah 
the azimuthal orders for degrees between 8 and 50; however, the coincidence of 
the I = 11 bump seems to have been merely a coincidence of noise, as it was 
not reproduced i n the early observations from the Big Bear Solar Observatory 



( LibbrechtlflQSl 



5.4 The SCLERA modes 



Hill et al\ (|l982i ) derived splittings from the SCLERA low-frequency peaks. 



and from those inferred a core rotating at 6 times the surface rate; however. 



WoodardI (jlQSJ) used ACRIM data to place an upper limit of 2.2 times the sur- 
face rate o n the interior rotation rate, inconsistent with these splittings. Later, 
HiH (ligss ") identified low-degree rota tional split tings in the five-minute band o f 



the SCLERA acoustic spectrum, but Libbrecht (|1986i) and Brown et al. i 19891 ) 



found that these results were inconsistent with the other evidence and were 
probably the result of misidentification of the modes. Given the complexity of 
the spectrum in question, whose derivation from measurements sampled at a 
few points on the solar limb made it difficult to separate out spectra of different 
degree, this seems a likely explanation. 

5.5 Low-degree acoustic mode splittings 1988 — 2002 

The next several years were active ones for low-degree helioseismology, with 
the development of the BiSON (Birmingham-based) and IRIS (based in Nice) 
networks. Together with the IPHIR instrument that rode the PHOBOS space- 
craft on its cruise phase to Mars, and the ground tests of the LOl (Luminosity 
Oscillations Imager) instrument that would later be mounted on the SOHO 
spacecraft, these brought a succession of estimates of the low-degree splitting, 
as summarized in Table [T] and Figure [141 In addition to the MDI instrument 
for medium and high-degree observations, the SOHO spacecraft carried both 
LOI and GOLF (Global Oscillations at Low Frequencies) specifically for ob- 
serving low-degree modes. Even though GOLF malfunctioned and could not 
be operated in its intended differential mode, instead being confined to making 
Doppler observations on one side of an absorption line, it provided some of the 
best available long-term, low-degree observations. 

The reported results show considerable variation, but apart from the early 
Tenerife result, which was based on much shorter and lower-duty-cycle observa- 
tions than most of the others, they all cluster around the surface rotation rate, 
some (particularly the IRIS results) pointing to a core rotation faster than the 
surface rate and some (in particular the BiSON results) to one substantially 
below it, perhaps as low as zero. As we approach the present time and the 
observation and analysis improve, the values tend to converge on a splitting 
quite close to that which would correspon d to the surface rate. Early in this 
period, there was room to speculate (e.g., Chaplin et al . 1996a[ ) that the dif- 



ferences reflected a temporal variation, but this could not explain away all the 
discrepancies. 
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Figure 14: I — 1 splitting estimates as a function of publication date. 



Table 1: Summary of ^ = 1 splitting measurements, 1988-2002 



Reference 


Project 


^I>'(/i,IIz) 


Comment 


Palle et al. (1988) 


Tcncrifc 


0.75 


summers of 1981-1986 


Toutain and Frohlich fl992) 


IPHIR 


0.563 ±0.017 


Intensity measurements 








on PHOBOS spacecraft 


Loudaeh et al. (1993) 


IRIS 


0.494 


Based on 3 low-frequency 








modes. 


Jimenez et al. (1994") 


Tenerife 


0.4768 ±0.0097 


Solar maximum 


Jimenez et al. (1994) 


Tenerife 


0.525 ±0.0127 


Solar minimum 


Toutain and Kosovichev (1994) 


IPHIR 


0.452 ±0.020 




Elsworth et al. ( 1995) 


BiSON 


0.42 ±0.02 




AoDOurchaux et al. (1995) 


LOI 


0.402 ±0.031 


I = 2 


ChaDlin et al. ( 1996a) 


BiSON 


0.415 ±0.006 




Lazrek et al. (1996) 


IRIS 


0.456 ± 0.012 




Lazrek et al. ( 1997) 


GOLF 


0.452 ± 0.014 




Gizon et al. ( 1997) 


IRIS 


0.456 ±0.010 




Bertello et al. (2000) 


GOLF 


0.436 ±0.009 




Bertello et al. (20001 


MDI 


0.447 ±0.011 


Asymmetric profile 


ChaDlin et al. (2001) 


BiSON 


0.435 ±0.0036 




Gellv et al. (2002) 


GOLF 


0.433 ±0.002 
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Figure 15: Power spectrum from 10 years of BiSON data, 1992-2002; the 
insets show the low-frequency end of the five-minute band (blue) and a single, 
rotationally spht I = 1 peak (red). 
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5.6 Pitfalls of low-degree splitting measurements 

Unfortunately, all the measurements described in Section [5.51 suffer from similar 
problems, as summarized below. 

1. The two components of the I = 1 mode are so close together (probably less 
than one microhertz apart) that they are resolved only for modes below 
about 2.2 mHz. This has implications for the measurements: 

(a) Estimates of the spli ttings of unresolved compon ents are highly prone 



to systematic error (jAppourchaux et aLl . l2000al ). 



(b) The components that can be resolved have small amplitudes (Fig- 
urell5p and therefore require both observations over extended periods 
and high signal-to-noise ratios. 

(c) On the other hand, these low-frequency modes have the advantage 
that they show very little frequency shift with the solar cycle, which 
simplifies the analysis of long time series. 

Even though the low-degree modes penetrate deep into the solar interior, 
they spend most of their time in the outer layers of the Sun and are not 
very sensitive to the core; conversely, estimates of the core rotation are 
very sensitive to small errors in the splitting measurements. 

In order to properly estimate the rotation profile in the deep interior it is 
necessary to combine the low-degree splittings with medium-degree ones 
in an inversion. However, because the low-degree modes are so few - 
a few dozen at most, compared to a couple of thousand medium-degree 
multiplets with tens of thousands of individual frequencies or coefficients - 
the need for extremely precise measurements is even more pressing. Also, 
combining data from different instruments with different systematic errors 
may cause problems, particularly if the observations were made at different 
epochs of the solar cycle. 



PointUabo ve was noted bylLoudaeh et aL J 1993) and Els worth et al. [ 19951 ). 
[ point m by iLoudagh et al\ (|l993l) and lLazrek et al\ (|l996h . who point out 



and 

that "An accuracy of about 30 nHz, or (1 year)"-'^ on the measurement of the 
I = 1 rotation splitting does not really permit, then, to discriminate between a 
solar core rotating twice as fast as t he rest or not ro tating at all!" An approach 
to addressing point [3] was made by iTomczvk et al. (fl995i) with the newly-built 



LOWL instrument, an imaging instrument optimized for lower degrees. They 
obtained splittings for 1 < / < 100, and inferred a rotational profile down to 
0.2 i?Q, finding a rotation rate that barely varied with radius between 0.2 i?0 
an d Q.QRp), apart from a low-s i gnific ance bump around 0.4 



I U.b Jtf?^, apart irom a low-s i gmhc ance bump around U.4it0. 
Eff-Darwich and Ko rzennik ('1998) further addressed point [3] when they com- 



bined results from several different instruments, including GONG, BiSON, MDI, 
and GOLF. They give a nice illustration of the tendency of higher-frequency 
low-degree mode splittings to be biased upward by the mode width, a point 
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that was further illustrated by IChaplin et al. ( 2001 ). and conclude that with 
the then-available data it is not possible to rule out fast rotation in the core 

bel ow 0.18^(71. 

ICharbonneau et al. I (|l998h used a genetic forward-modeling approach to an- 
alyze the LOWL data, with results favoring a rigidly-rotating core. 



5.7 A new millennium for low-degree helioseismology 

Starting around the turn of the century, there was a move towards more collab- 
orations and comparisons between different projects in an effort to understand 
the systematic errors and better constrain the solar core dynamics. By this time, 
multi-year observations were available from GONG and the SOHO instruments, 
as well as good-quality observations from BiSON stretching back to 1991. 

Chaphn et M (|l999l) combined the LOWL higher-degree splittings with the 



very precise low-frequency BiSON splittings for the lowest-degree modes, and 
concluded that the data were consistent either with rigid rotation or with a 
slight downturn in the rotatio n rate in the c o re (the latter being at best a 1- 
cr result); on the other hand, Corbard et al. ( 1998bh had used a very similar 
analysis of GOLF and MDI data to deduce a slight increase in the rotation 
rate below 0.25 i?©, but Garcia et al. ( 2003[ ). also using MDI and GOLF data, 
obtained rather low splitting values from a 2243-day time series and tentatively 
co ncluded that they c ould rule out a high rotation rate i n the core. 



^ Eff-Darwich et ali (|2002,) , following on from the work of lEff-Darwich and Korzennik 
(ll998|), again combined BiSON, GOLF, GO NG and MDI data and found a very 



small downturn in rotation in the core, while ICouvidat et ali (|2003l ) found a flat 
rot ation profile down to 0.2 Rq using combined GOLF, MDI and LOWL data. 
Fletcher et all (|2003h investigated the problem of fitting the poorly-resolved 



higher-frequency low-degree mode splittings to integrated-sunlight observations 
such as those from BiSON. Using genetic fitting algorithms, they were able 
to reduce, though not eliminate, the bias towards higher splittings for these 
fits. They also found, in common with previous work, a strong anticorrelation 
between the estimated splitting value and its formal error, which would tend to 
ca use overestima t ed sp littings to be more heavily weighted in inversions. 

Garcia et al. (|2004l) considered two years of "sun-as-a-star" observations 



from early in the solar cycle, obtained from GOLF, GONG, MDI, VIRGO and 
BiSON, and were able to extract not only sectoral splittings but also and 
coefficients from the data, suggesting that it may be possible to infer differential 
rot ation even in s t ars fro m which we will never have resolved data. 

IChaplin et al. (2004) used artificial data to address the question the de- 
tectability of a rotation-rate gradient in the core. They concluded that, based 
on the best available data from ten years of observations, the difference between 
the rotation rate at 0.1 Rq and 0.35 Rq would be detectable only if it exceeded 
11 nHz. 

Chaplin et al. ( 20061 ) carried out an exhaustive "hare-and-hounds" exercise, 



in which one participant (the "hare" supplies the same set of artificial data to 
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the others, the "hounds," who then apply their various fitting methods with- 
out knowing the "true" answer, and compare the results. They obtained good 
agreement between the different techniques for I — 1, but systematic differences 
for the / > 1 splittings, which are attributed to different assumptions about the 
relative heights and spacing of the non-sectoral (|m| < /) components. 



5.8 Summary of the acoustic-mode results 

To summarize, the best evidence we have so far seems to imply that the ro- 
tation rate between about 0.2 i?Q and the base of the convection zone is most 
likely approximately constant with radius and spherically symmetric. It is not 
possible to rule out a different rotation rate for the inner core, but there is no 
evidence from p-mode observations to support such a difference. Between about 
0.2 i?0 and the base of the tachocline, n o significant departure fro m rigid-body 
rotation has been found. As discussed by EfF-Darwich et al. ( 20021 ). for example 



the available constraints already seem to rule out the simplest models of hydro- 
dynamic spin-down, which would show a detectable increase in the rotation rate 
below 0.3 Rq. Understanding both of the relationship between p mode splittings 
and the interior rotation, and of the care needed to measure them, has greatly 
advanced since the early days of helioseismology, but the rotation rate of the 
innermost nuclear-burning core remains uncertain. 



5.9 Gravity modes 



One possible way to improve the constraints on the core rotation would be to 
use g modes, or gravity waves, instead of p modes. Because these modes have 
their greatest amplitude in the solar interior, they should be much more sensi- 
tive to the core properties. Unfortunately, they also have very small amplitudes 
at the surface. The history of helioseismology is littered with unconfirmed re - 
ports of ff-mode ide n tification; see, for exa mple, Delache and Scherrei ( 19831): 
Van der Raavl(|l988f) : [Thomson et a/.l (|l995l) . and the review bylHilLeTa/] \\MVd) 
The most promising recent work has been carried out using long ti me series from 
the GOLF instrument aboard SOHO. Appourchaux et al. ( 2000bf ) placed an up- 
per limit of lOmm / s on q -mode amplitudes based on two years of observations, 
and .Gabriel et al\ ( 2002) reduced thi s limit further, to 6mm/s, using 5 years 
of data. Most recently, iGarcia et !all (|2007l ) report finding a pattern of peaks 
with constant spacing in period corresponding to the model-predicted spacing 
for I = 2 g modes with SI = 0, (5n = 1, and with a splitting that they interpret 
as corresponding to a core rotation rate of 3 - 5 times the surface rate; however, 
this is still a prelimina ry result in ne e d of c onfirmation. 

In a related paper, iMathur et al\ (|2007( ) point out that the current predic- 
tions for low-order g-mode frequencies are much more consistent than was the 
case a decade earlier, resulting in a period for the fundamental g-mode between 
34 - 35 minutes. This finding does make one wonder about the usefulness of 
the g-mode observations for discriminating among models; on the other hand, 
it lends somewhat more credence to the current identification. 
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6 The Tachocline 



While the bulk of radiative interior appears to rotate almost as a solid body, 
the base of the convection zone at 0.71 Rq coincides with a region of strong 
radial shear, above which the convection zone exhibits a differential rotation 
pattern that depends strongly on latitude and only weakly on depth. This 
shear la yer is known as the tachocline, a term introduced to the literature by 
ISoiegel and ZahnI (Il992h . who att r ibute to D.O. Gough the correction of the 
earlier term "tachycline" ( Spiegel . Il972l ). As is evident from the date of the 
latter reference, the notion of a shear layer at the bottom of the convection zone 
had been present in models for some time prior to its observational discovery, 
though its exact location was somewhat uncertain. 

The existence of a layer of radial shear around the base of the convection 
zo ne, with appr o ximat ely solid-body rota tion below it, was fi r st dem onstrated 
bv lBrown et al. ( 19891 ). using the data of Brown and Morrow ( 1987 ): however, 
the significance of their results was quite low and the y were at pain s to po i nt out 
that other interpretations of the data were possible. iDziembowski et ^ (Il989l ) 
used BBSO data to improve the picture of rotation at the base of the convec- 
tion zone, again finding that the low-latitude rotation rate increased, and the 
high-latitude rate decreased, towards a common value at the base of the convec- 
tion zone. The position of the bas e of the convection zone was determined by 
IChristensen-Dalsgaard et a l. ('1991") using so und-speed inversions of helioseismic 
freque ncies from the work of D uvall Jr et al\ (il988) andLibbrecht and KaufmanI 
( 1988f ): their value of 0.713 i?©, confirmed by Basu and Antia (1992), has been 
accepted ever since. 

The discovery of this shear layer (as pointed out by Brown et al.) olfered a 
solution to the puzzle of the apparent absence of a radial gradient of rotation 
in the convection zone that could drive a solar dynamo, leading to speculation 
that the dynamo must operate in the tachocline region instead of in the bulk of 
the convection zone. 

The tachocline lies in the region where modes of Z « 20 have their lower 
turning points, and the resolution of the inversions is quite low - about 5-10% 
of the solar radius in the radial direction. The thickness of the shear layer 
is therefore likely not to be resolved in inversions, and some ingenuity (and 
forward modeling) is required to estimate it and account for the effect of the 
finite- width averaging kernels in smoothing out the inversion inferences. The 
results of various efforts to parameterize the tachocline shape at the equator are 
summarized in Tabic [21 They mostly concur in placing the centroid of the shear 
layer slightly below the seismically-determined base of the convection zone, and 
its thickness at arou nd 0.05 i?Q. The largest value for the thickness, that of 
IWilson et all (Il996bl) . was obtained using forward calculation and direct com- 
bination of_S£litting_coefficien rather than a true inversion, while the very low 
value of Corbard et al. ( 199 9^1 was obtained using an inversion technique specifi- 
cally designed t o allow a discontinuous ste p in the rotation rate at the tachocline. 
The analysis of lElhott and Goughl (|l999l ) was somewhat different from the oth- 
ers, in that it involved calibrating a particular model of the tachocline against 
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Table 2: Tachocline radius r and width F. 



Reference 




(Tr/ Rq 


T/Rq 


crr/R@ 


Project 


Kosovichcv (1996) 


0.692 


0.005 


0.09 


0.04 


BBSO 


Wilson et al. fl996a) 


0.68 


0.01 


0.12 




BBSO 


Basu (1997) 


0.705 


0.0027 


0.0480 


0.0127 


GONG 


Antia et al. i 1998) 


0.6947 


0.0035 


0.033 


0.0069 


GONG 


Corbard et al. ('1998a) 


0.695 


0.005 


0.05 


0.03 


LOWL 


Corbard et al. (1999) 


0.691 


0.004 


0.01 


0.03 


LOWL 


Charbonneau et al. (1999) 


0.693 


0.002 


0.039 


0.002 


LOWL 


Elliott and Gough fl999) 


0.697 


0.002 


0.019 


0.001 


MDI 


Basu and Antia (2003) 


0.6916 


0.0019 


0.0162 


0.0032 


MDI, GONG 



th e inferred sound-sp eed r ather than agai n st a r otation profile. 

Antia et al\ (|l998l ) and ICorbard et ail (|l999l ) found no significant evidence 



for a variation in the posi tion or thickness of the tachocline with latitude, but 
Charbonneau et al. I (11999') found a significant prolateness, wit h the tachocline 
(0.024 zfc0.004) Rq shallower at latitude 60° than at the equator. iBasu and Antia 
(l2003h also found a slightly thicker and shallower tachocline at high latitudes, 
and speculated that the tachocline location might be discontinuous at the lati- 
tude (around 30°) where the shear vanishes and changes sign. 



6.1 Models and the tachocline 

Even the most generous estimates for the observed tachocline thickness are small 
enough to pose an interesting theoretical question: what prevents the shear from 
spreading further into the radiative interior, destroying the observed uniform 
rotation? The literature on tachocline modeling is extensive, far beyond the 
scope of this review. In brief, three main c andidate mechanisms have been pro- 
posed: turbulent flows ([ S piegel and Zahnl . Il992) : "fossil" magnetic fields (e.g., 
Gough and Mcintvre 19981 ): an d gravity waves, known to observational helio- 
seismologists as g modes (e.g., Zahn et al. 1993), but all these scenarios have 
problems when considered as the sole mechanism. Recent advances in comput- 
ing have made possible detailed three-dimensional simulations to explore these 
issues, but these models have not yet been able to reproduce a self-sustaining 
tachocline. For a review from a modeler's perspective, see iMiesch (2005). Also, 
a variety of d i scussi ons of tachocline models are collected in the book edited by 
Hughes all (|2007l ). 
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7 Rotation in the Bulk of the Convection Zone 



The surface differential rotation, with the equator rotating faster than the poles, 
was known from, for example, sunspot tracking, long before helioseismology 
opened up the solar interior. Most models in the pre-helioseismology era pre- 
dicted or assumed a rotation rate constant on cylinders parallel to the axis of 
rotation. This is a consequence of the so-called Taylor-Proudman constraint, a 

we ll-known result in fluid dynamics. 

Duvall Jr and HarvevI (|l984l) ; lDuvall Jr et all (|l984l ) observed from the South 



Pole, using only sectoral modes; their instrument used intensity images in a Cal- 
cium absorption line, scanning rather than imaging the whole Sun at once. Their 
main conclusion was that, "Most of the Sun's volume rotates at a rate close to 
that of the surface." 



Brownl (|l985f ) had a different instrum ent, th e Fourier Tachometer, which 



produced 100 x 100 pixel velocity images. iBrownf s initial crude analysis of five 
days of data used cross-correlation, and expanded the multiplet frequencies using 
low-order polynomial fits; the results showed little sign of any depth variation 
in the differential rotation. 

), again using data from South Pole observations but 



Duvall Jr et al. 



now covering the full range of azimuthal orders, found values of the coeffi- 
cient (the first-order measure of differential rotation) consistent with the surfac e 
rot ation and rather larg e r than was consistent with the results of lBrownlljlQS^h . 

iBrown and Morr"ow ( 1987 ). with 15 days (not all consecutive) of Fourier 
Tachometer data, could not distinguish between rotation on cylinders and latitude- 
dependence, but found that there was definitely less differential rotation in the 
radiative in terior below the conve ction zone; their 03 values were now closer 
to thos e of Duvall Jr et al. ( 1986h . and they declared the previous ones erro- 



Brown et al. carried out a much more detailed analysis of the 

Brown and Morrow (198Z) data, strengthening the evidence for mostly depth- 



neous. 



independent rotation in the convection zone, as shown in Figure [TBI 

Both the South Pole observations and those of Brown and collaborators 
were relatively noisy and of poor resolution; although they strongly hinted at 
a picture with little radial differential rotation in the convection zone and little 
differential rotation at all below it, other interpretations were possible. 

iLibbrc cht (1989,) published splittings from 100 d ays of BBSO observ ations 
in summer 1986, broadly confirmin g the results of Brown et al. ( 1989f ) with 
substantially smaller uncertainties. iDziembowski et all (|l989r ) inverted these 
data, and inferred a sharp radial gradient at the base of the convection zone 
and roughly constant rotation at each latitude above that. They also found a 
bump in the rotation rate in the middle of the convection zone, to which we 
will return below. Other inversions of the s ame data s et were presented by 
Christensen-Dalsgaard and Schou ( 1988h and iLibbrechtj (1988b), with s i milar 



resul ts, though not all the early inversions (c.f. iKorzen nik et al. Il988t ISekiil 



19911 ) produce such recognizable results; this may be an example of the difficulty 



of using OLA-type techniques for data with insufficient higher-degree modes. 
Another (2D OLA) inversion of these data, shown in Figure flTl was carried out 
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Figure 16: Summary of rotation inferences from lBrown et al. ( 1989I ). reproduced 
by permission of the AAS. 



bv lSchou e^"aLl(|l992h . who illustrated their averaging kernels; these were rather 
broad, but adequate to rule out a rotation-on-cylinders model. This paper was 
also the first to make the important point that the so-called "polar" rotation 
rat e inferred from inver sions is actually localized somewhat away from the pole. 
Gough et al. I (|l993l) continued to challenge the observers to completely ex- 



clude rotation on cylinders, pointing out that it was possible to construct a cylin- 
drical m odel that satisfied the constraint of the BBSO data, but lSchou and Brown 
(19943) showed that such a model could not be made consistent with both the 
Fo urier Tachometer data and the gravitational stability of the rotating Sun. 

Bachmann et al. ( 19931 ) analyzed Fourier Tachometer observations from 1989 
and pointed out a "wiggle" in the splitting coefficients at z^/L ~ 40/iHz, (corre- 
sponding to a turning-point radius of about 0.85 -/?©); attributed to daily modu- 
lation of the observations, this now well-known effect accounts for the "feature" 
seen in the middle of the convection zone in many inversions of single-site data. 

Better data, with long time series free from daily modulation, were obvi- 
ously needed before much more progress could be made, and with the advent 
of the GONG network in 1995 and the MDI instrument aboard SOHO in 1996 
such data became avai lable. Preliminary r otation profiles w e re pr esented by 
Thompson elo/] (|l996l ) for GONG and by iKosovichev et al\ (|l997l ) for MDI, 
both showing the now familiar pattern of almost-constant rotation in the con- 
vection zone, with shear layers both at the base of the convection zone and 
below the surface. 



Schou et al\ (|l998f ) carried out a comprehensive analysis of the rotation pro- 



file based on the first 144 days of observations from MDI, using and comparing 
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several different rotation inversion techniques with an input data set consisting 
of coefficients up to 035 for p modes up to ^ = 194 and / modes up to / = 250. 
They were able to obtain consistent and robust results from the surface to about 
0.5 i?0 at low latitudes; at higher latitudes the domain of reliability was shal- 
lower. Roughly speaking, the inversions could not be well localized within about 
0.2 i?Q of the rotation axis. The results (Fig. [T5)l showed that the rotation in 
the bulk of the convection zone, below 0.95 i?o, had a slow increase with radius 
at most latitudes, but was definitely incompatible with rotation on cylinders. 




Figure 18: Four inferre d rotation profiles from the first 144 days of MDI ob- 
servations (ISchou e< a/.l,[l99 8); (a) 2DRLS, (b) 2D0LA, (c) IDxlD SOLA, (d) 



1.5d RLS, from Schou et al\ (| 19981 ) . reproduced by permission of the A AS. 



7.1 The "polar jet" 



In addition to the other general features described here. lSchou et aLl found some 
evidence for a "jet" of faster rotation at about 75° latitude and 0.95 i?©; al- 
though this was more obvious in some inversions than in others, it did seem 
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to have a signature in the coefficients themselves (see also Howe et aZlll998f ). 
However, this feat ure was not reproduced in inversions of GONG data (e.g., 
Howe et oL ^ 2000b), or even in inversions of MDI data analyzed with the GONG 
pipeline (jSchou et aLl . 120021 ). and it is now believed to be an artifact related to 
the MDI data analysis. 



7.2 GONG/MDI comparison 

Once both GONG and MDI had been running for a few years, it became ev- 
ident that the two projects were producing inferences of the interior rotation 
profile that were different in some si gnificant details, p articularly at high lat- 



itudes within the convection zone. ISchou et al\ (j2002r ) carried out a careful 



comparison, taking data from three epochs at different phases of the solar cy- 
cle from each project and deriving rotational splittings or splitting coefficients 
from each, both with the usual algorithms and with those regularly used for the 
other project's data, before using both RLS and OLA inversions. The results 
clearly showed that most of the discrepancies arose from the analysis pipelines 
rather from the data themselves. The "CA" peak-fitting algorithm used for the 
MDI data was able to extract modes from the GONG data to somewhat higher 
degrees and lower frequencies than the "AZ" algorithm could manage with ei- 
ther GONG or MDI input data. However, for both MDI and GONG data, the 
"CA" algorithm introduced an anomaly in the splitting coefficients centered at 
around 3.3 mHz, which in turn caused the inversion inferences to show a higher 
rotation rate deep in the convection zone at higher latitudes. Excluding these 
data brought the GONG and MDI data (analyzed with the "AZ" and "CA" 
pipelines respectively) into much better agreement, at the cost of somewhat 
degraded resolution. Restricting both data sets to the common mode set below 
3 mHz reduced the discrepancies even farther, but did not remove the "jet" in 
the MDI data. Since the "jet" feature was only seen in the MDI data analyzed 
with the CA pipeline, however, the authors concluded that this feature was 
probably spurious. 



7.3 Slanted contours 

Although much of the debate in the early 1990s centered on discriminating be- 
tween rotation constant on cylinders and rotation constant along radial lines 



Oilman and Howe 



neithe r picture gave a co mplete description of the data, 

( 2003h : lHowe et al\ (120051 ) jointed out that the differential rotation in the bulk 
of the convection zone, at least at low- to mid-latitudes, could be quite well 
described by saying that the contours of constant rotation lay at about a 25° 
angle to the rotation axis, as illustrated in Figure [T91 

Figure [201 compares idealized rotation profiles for the cylindrical, radial, and 
slanted-contour configurations. 
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Figure 19: Mean rotation profile from GONG data; contours of constant ro- 
tation (left) , show ing lines at 25° to the rotation axis as dashed lines, after 



Howe et al\ (120051), and c uts at constant latitude as a function of radius (right) , 



after iHowe et aLr(|2000bf) 



7.4 Polar rotation 



Another interesting feature revealed by the early GONG and MDI observations 
( Schou et ai . 1998 : Birch and Kosovichevl Il998l ) was that, while the surface 



rotation rate was mostly well described by the usual three-term ex pansion in 
the c osine of the colatitude 9, i}{d) — A + B cos^ 9 + C cos^ 9, (e.g., ISnodgrasa 



19841 ) the rotation rate close to the poles was significantly slower than that. 
The authors speculated that this might be a result of drag from the solar wind, 
and that the effect might therefore disappear or become less marked at epochs 
of higher activity. In fact, though the inferred high -latitude rat e did s peed up 
during solar maximum - as seen, for example, in iHowe et al. (|2005l) and in 
Figure US] - it remained at all times lower than the extrapolation of the three- 
term fit. 



7.5 Models and rotation in the convection zone 

The interior rotation is only one part of the complex system that drives the solar 
cycle, but it is perhaps still the easiest part to measure in the solar interior; the 
meridional circulation can be directly measured only in the shallower subsurface 
layers, and buried magnetic fields can at best only be inferred indirectly. The 
differential rotation in the convection zone must arise from the interaction of 
convection cells and Coriolis forces, with the meridional motions playing an 

important part. 

Early depictions of the solar dynamo (see, for example, Kohleil[l974 : Durnev 



m required a rotation rate increasing inward, and a meridional ^o^ rl^t 



the poles and sinking at the equator, in order to drive the solar cycle migration 
of the activity belts in the observed sense. This picture, taken together with 
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Figure 20: Idealized rotation profiles for rotation constant on cylinders (left), 
radial lines (middle) and lines at 25° to the rotation axis (right). The top 
row shows contours of constant rotation, while the lower row shows rotation 
rate as a function of radius at constant latitude for latitudes at 15° intervals 
from the equator (top) to 75° (bottom). The rotation rate is matched to the 
GONG inferences at 0.99 Rq and smoothed to simulate the broadening effect of 
inversion resolution on the tachocline; the near-surface shear was not included. 
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rotation on cylinders, would have meant that the observed surface differential 
rotation was a superficial phenomenon, with the dynamo operating in the un- 
observable deeper layers. At this stage, there does not seem to have been a 
clear distinction made between the direction of the meridional circulation at 
the surface and the direction of migration of the magnetic activity belts dur- 
ing the solar cycle, which are of course now understood to operate in opposite 
directions; the poleward meridional flow at the surface was first measured by 



Duvall Ji (197' 



The models of Glatzmaiei (|l985l) and lGilman and Milled (|l986[ ). which were 
among the first numerical simulat i ons o f solar rotation and the dynamo, have 
been cited, for example bv lWilsonI (|l992l ) as dating from "Prior to the advent of 
helioseismology," but this is not q uite correct. In fact , both these papers refer 
to the Duvall a nd Harvey dat a, and Gilman and Miller ( 1986f ) also mentions the 
observations of Brown ( 1985[ ). suggesting that the model results could be consis- 
tent with the helioseismic observations if there were a layer of inward-increasing 
velocity below the surface and above the domain of the simulation. The simula- 
tions in both case s, like their precursors ove r the previous several years such as 
that described bv iGilman and Miller (1981), produced rotation approximately 
constant on cylinders and increasing outward, which would result in a dynamo 
wave propagating poleward if the dynamo were operating in the bulk of the 
convection zone. The main message that modelers in the late 1980s seem to 
have taken from the observations was that th e rotation rate was i ncreas ing out- 
ward, in agreement with the simulations of Gilman and Miller ( 19861 ) but in 
disagreement with the a-effec t dyna mo picture, which required a rotation rate 
increasing inward; see iParkeri (jl987l ) for a review re presenting a theorist ' s per - 
spective on the state of play at this stage. This led Gilman and Milleil (1986 *) 
to su ggest (not for the first time; see also, for example, I Galloway and Weisa 



198lh that the dynamo might be operating in a thin layer at the bottom of the 



convection zone; this speculation was further reinforced by the later helioseismic 
inferences that clearly showed this shear layer, or tachocline (see Section[S]) and 
the approximately radial configuration of the rotation in the convection zone. 

Even quite recent global simulations of convection ( Brun et al\ . 2004 . for ex- 
ample, ) still show some tendency towards ro tation on cylinders, but the higher- 
resolution calculation of iMiesch et ali (j2008r ) mostly eliminates the cylindrical 
effect and produces a rotation pattern, based on giant convection cells, that after 
suitable temporal averaging looks quite solar-like, as illustrated in Figure [21] 
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Figure 21: Three tem porally ave r aged rotati on profiles from the spherical- 
shell simulations of (a) Brun et al\ ( 2004 ). (b) iBrowning et al. ( 20061 ). and (c) 
( 20081 ). reproduced by permission of the AAS. 



Miesch et al. 
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8 The Near-Surface Shear 



One persistent puzzle in the measurements of rotation at the photosphere had 
been that direct Doppler measurements consistently gave somewhat slower ro- 
tation rates than th e meas urements made by tracing s urface features. For ex- 
ample, iBrown et summarized the results of ISnodgrasd (|l983l \l984) 

as 

462 - 74^2 „ 53^/nHz (25) 



2tt 



for magnetic features and 



^ = 452 - 49Ai^ - 84^^nHz 

ZTT 



(26) 



for the surface plasma, respectively, wh ere ii i s the sine of the latitude. For 
an overview of such measurements, see iBeck (l2000l) . The usual explanation 
for the discrepancy is that while the Doppler techniques measure the velocity 
at the surface, the tracers such as sunspots are ancho red in a faster-rotating 
layer deeper down. For example, Gilman and Foukail |l979) noted that the 
observations implied a subsurface shear layer and suggested that this might 
arise from angular momentum conservation in the supergranular layer. 

An extremely early at tempt to measure the subsurface rotation was made 
bv [Rhodes Jr et all (|l979f) . when the identification of the 5-minute oscillations 
with p modes was still a relatively recent discovery. These authors used high- 
degree modes, probing about the upper 20 Mm (0.03 i?©) of the convection 
zone, and detected an inwards-increasing gradient. If these measurements are 
reliable, they represent the first detection of the subsurface shear. However, 
most of the early helioseismic measurements of the internal rotation profile 
were restricted to a degre e range that did not al low the near-surface shear to 
be resolved in inversions. Rhodes Jr et al. ( 1990h . attempting to measure the 
rotation in the bulk of the convection zone, also saw hints of a gradient, opposite 
to tha t seen at the base of the convection zone, below the surface, and I Wilson 



1992) used forward c a lculat ion techniques on the data of iBrown and Morrow 



mS) and iLibbrechtJ (Il989t l to deduce that the rotation rate must increase 



inward immediately below the surface. We should remember, however, that at 
this time the picture of the internal rotation profile was not as clear as it is 
today, and it is not always obvious whether interpretation of the observations 
as gradients of rotation refers to the near-surface shear, the shear at the base of 
the convection zone, or some unresolved amalgamation of the two. Wilson, for 
example, was not arguing for a near-surface shear layer but against the model 
with rotation constant along radii. 

With the advent of GONG and MDI, measuring modes to higher degrees 
than had previously been possible, the near-surface she a,r could be seen in globa l 
inversions; it is visi ble in the early re s ults p resented bv [Thompson et~al. i 19961 ) 
for GONG and by iKosovichev et all (|l997[ ) for MDI, in both cases apparently 
changing sign at higher latitudes. 
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Schou et al\ (|l998f ) found clear evidence of the near-surface shear in inver- 
sions of MDI data. All the inversion methods agreed well on the shear at low 
latitudes, but at high latitudes the picture was complicated by the proximity of 
the submerged "jet" feature and the methods agreed less well. The disagreement 
may have been partly due to systematic errors in the splitting coeffic ients . In the 



compa risons of MDI and GONG data and analysis carried out by ISchou et~ai. 



(2003), the high-latitude reversal of the shear is seen only in data analyzed with 
the "CA" pipehne; this may be partly because the "AZ" pipeline mostly fails 
to recover the splittings of the (narrow, low-amplitude) /-mode peaks, but the 
reversal persists in the MDI data even for the restricted common mode set. 
The near-surface shear (do wn to about 15Mm) was studied in detail by 



Corbard and ThompsonI (|2002[ ). using / modes from MDI data. They measured 
the slope of the rotation rate, close to the surface at low latitudes, as about 
—400 nHz/i?0, decreasing to a very small value by about 30° latitude and 
possibly reversing in sign at higher latitudes (though this result, seen in only 
the outer 5 Mm, was dependent on only the highest-degree modes, those with 
/ > 250). The low-latitude rotation rate was found to vary almost linearly with 
depth in the subsurface region, while if angular momentum was conserved in 
parcels of fluid moving with respect to the rotation axis, it would be expected 
to vary with the inverse square of the distance from the axis. 

The near-surface shear is also accessible to the methods of loc al helioseis- 
mology, at least for latitudes below 50-60°. Basu et all ( 19991 ) and Howe et al 



(|2006ah compared results from local ring-diagram analysis and global inversions 
and found, at latitudes < 30°, quite good agreement between the dV,/dr values 
obtained from local and global inversion results. However, although the slope 
from local measurements does show some variation with latitude (Figure [22|). it 
by no means vanishes at 52.5°, the highest latitude at which the measurement is 
made. The ring-diagram re sults allow us to c onsider the northern and southern 
hemispheres separately, but iBasu et ali ([19921) found very little difference in the 
shear between the two hemispheres. 

Some attempts have been made to use the near-s urface shear to drive or at 
least contribu t e to a solar dynamo, for example by iBrandenbure (|2005[ ). but 



Dikpati et ali (|2002i ) showed that any dynamo contribution from the shear of 



the outer layers could only provide a fraction of the effect needed to power the 
solar cycle. 
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Figure 22: Radial variation of the mean rotation rate after subtraction of the 
tracking rate, for global inversions (blue) and north -south averaged local in- 
versions of MDI (green) a nd GONG (red) dat a at latitudes 0° (a), 15° (b), 30° 
(c) and 45° (d); similar to lHowe et M (|2006al ). 
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9 The Torsional Oscillation 



The so-called Torsional Oscillation is a pattern of migrating bands of faster- and 
slower-than-average zonal (i.e., parallel to the equator) flow associated with the 
equatorwar d drift of the activity belts during the solar cycle. It was first de- 
scribed by [Howard and Labonte (|l98d) . who used twelve years (1966-1978) 
of full-disk velocity observations from the 150-foot tower at the Mount Wil- 
son observatory and found evidence of a pattern of flow bands migrating to- 
wards the equator; the greatest concentration of active regions is associated 
with the poleward edge of the main equatorward-moving band. They initially 
interpreted the high-latitude variations as consisting of bands of faster rotation 
starting at the poles and tak ii ig a full 22-year Hale cycle to drift to the equator. 
Scherrer and Wilcoxl (|l980a[) : IScherrer et al\ (|l980l ). observing at the Stanford 
Solar Observatory, found no evidence of changes in the equatorial rotation rate 
for data from 1976- 1979, but as this period was close to a solar minimum, and 
the resolution of the Stanford instr ument was not high, this is neither surprising 
nor in consistent with the results of Howard and Labonte . Labonte and Howard! 
(| 19821 ) note that IScherrer and Wilcoxl (|l980bl ) fat a AAS meeting), had "con- 
firmed the existence of the global velocity field," though this is not apparent 
from the latter's published abstract. 

A somewhat different pattern of velocity variations is seen when magnetic 
features rather than Doppler measureme nts are used t o determine the surface 
rotation rate, as described for example bv lKomm et "aLl (fl993 a) , who found that 
the pattern derived from magnetograms lay equatorward of that from Doppler 
measurements, with the slower-than-average bands coinciding with the zones of 
greater magnetic flux. 

Mount Wilson Doppler observations since 198 6, clearly showing th e pattern 
of migratin g zonal- flow bands, were presented by lUlrichI (Il998l . l200l!) : see also 
Howe et MI (l2006al) for updated results. The bands extend over about 10° in 
latitude, and have zonal velocities a few meters per second faster or slower than 
the surrounding material, corresponding to excess angular velocity of less than 
0.5% of the overall rotation, or a few nanohertz. 



9.1 Early helioseismic measurements 

The first hints of the signature of the migrating fiow bands in helioseismic data 
can be seen in the BBSO data (jWoodard and Libbrechtl . Il993l ). as was pointed 
out bv lHowe et al. ( 2000c ). but these measurements do not give mu ch informa- 
tion on the radial extent of the flows. Kosovichev and Schou! ([l993) found evi- 
dence of the flows, a few meters per second faster th an the general rota tion pro- 
file, in /-mode measurements from early MDI data; Giles et al. i 1998^) found a 
similar pattern using the time-distance techn i que of local helioseismology, while 
Schou and The SOI Internal Rotation Team! (|l998! ) and ISchoul (|l999[ ) clearly 
showed that these flows were migrating in a manner consistent with the Mount 
Wilson Doppler observati ons. The flrs t radial ly-resolved evidenc e of zonal flow 
migration was reported bv lHowe et oD (|2000c! ) for GONG and bv lToomre et al 
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(|2000l) for MDI, while iHowe et all (|2000al) combined MDI and GONG data and 



concluded that the equatorward-migrating part of the flow pattern (at latitudes 
below about 40°) pene trated to at least 0.92 i?© (56 Mm below the surface) . 
Antia and Basul (|2000D also reported similar findings. lAntia and Basul (|200l[ ) 



studied the evolution of the variations poleward of 50°, which had much higher 
amplitudes than the equatorward-moving flows and which showed signs of prop- 
agating poleward over time. The larger amplitude of the high-latitude signal 
may be related to the smaller angular momentum closer to the rotation axis. 

9.2 Recent results 

As more data accumulated, the signature of th e torsional oscillation p attern in 
the helioseismic observations became clearer. Vorontsov et ^ (|2002l) studied 



the evolution of the flows in MDI data from 1996 through 2001. They con- 
cluded that at least the high-latitude region of changing rotation involves the 
whole depth of the convection zone. The results on the radial extent of the flows 
at lower latitudes were less clear, with evidence that the bands of slower rota- 
tion might penetrate close to the base of the convection zone, while the bands 
of faster rotation appeared to reach about 0.9 Rq but no deeper. Another inter- 
esting feature of that paper was the introduction of the use of 11-year sinusoids 
to characterize the variation of the rotation rate at any given location. This 
innovation had the useful effect of clarifying the pattern, making obvious the 
poleward propagation of the high-latitude flows even with data from little more 
than half a cycle. The existence of a weak third-harmonic component to the 
ele ven-year cycle, h o wever , was not confirmed in later work. 

Basu and Antia (|2003h found similar results in MDI and GONG data up to 



2002, as seen in Figure [23l These results also hint at another subtlety; at low 
latitudes, the phase of the flow pattern is not constant along radial lines. In fact, 
the variation in the lower part of the convection zone appears to lead that close 
to the surface by a year or two, with the low-latitude band of faster rotation fol- 
lowing roughly th e same 25° slant as the rotation contours. This tendency was 
further studied by iHowe et al\ (|2005l . [2006b) , who compared inversions of MDI 



and GONG data with forward-modeled profiles based on different flow conflgu- 
rations, including some derived from dynamo models. Although some detail was 
lost and distorted due to the resolution and uncertainties in the inversions, the 
authors were able to conclude that the low-latitude branch probably penetrates 
through much of the convection zone, but is sufficiently displaced in phase at 
greater depths that the correlation between the surface pattern and that deeper 
down almost vanishes. In this work, the 11-year sinusoid analysis showed evi- 
de nce of a second-h armonic component rather than the third harmonic reported 
bv I Vorontsov et al. . 



Figures [25l [26l a nd [27l show the variations in rotation rate, based on the 



results and figures in iHowe et a,l\ (l2005l . l2006bh . but brought up to date with 



the most recent GONG and MDI observations available at the time of writing. 
The plots were prepared using the same 2-D RLS inversion codes for both MDI 
and GONG medium-degree data, and 2-D SOLA for MDI, that were used for 
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Figure 23: Contour diagrams of constant rotation velocity residuals at 0.98 i?©, 
o btained using two d imensional RLS inversion of the GONG data, from 



Basu and Antial (|2003f ). reproduced by permission of the AAS. 
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Figure 24: Zonal flow pattern derived from MDI /-mode measu rements, with 
smooth profile subtracted. Based on a figure from ISchoul (|l999l ). updated and 
used by kind permission of J. Schou (2008, private communication.) 



the work of lHowe et ali (j2000a[ l and the other related papers. Figure [28] shows 
the phase and amplitude profiles for 11-year sine functions fitted to the rotation 
variations. 



9.3 Local helioseismology and the torsional oscillation 

The torsional oscillation pattern, at least at lower latitudes and closer to the 
surface, is also suitable for measurements using the techniques of local helio- 
seismology, in which short- wavelength, short-lived waves are used to infer the 
structure and dynamics of localized areas of the Sun. Because these waves do 
not penetrate very far below the surface, such techniques are restricted to the 
outer few megameters of the solar envelope, but this region can be studied in 
much greater detail and with shorter averaging times than is possible with global 
hel ioseismol og y. 

Basu and"^A^rtial (|2000h dete cted the zo nal flow migration using MDI data 



and the ring-diagram technique (|Hi lil. ll988h . in which the displacement of three- 



dimensional acoustic power spectra derived from small areas of the solar disk is 
used to infer horizo ntal flows in both the zonal and meridional directions. Later, 
Haber et all (|2002h measured both the zonal flows and a correspo nding mod 



ulation of the meridional flow pattern, as seen in Figure [29] (left). iBeck et al 



(2003), using the time-distance technique, which considers the correlations be- 



tween oscillations at spatially separated locations, also found bands of merid- 



51 




a: h: 




■■1. ri:-. I A. '^1 

r ■ ^ ■ - ' r - J ^ 



Figure 25: Rotation rate after subtraction of a temporal mean at each loca- 
tion, as a function of latitude and time at selected depths, for OLA (top) and 
RLS (middle) inversions of MDI data, and for RLS inversions of GONG data 
(bottom). 
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Figure 26: Rotation rate after subtraction of a temporal mean at each loca- 
tion, as a function of depth and time at selected latitudes. Latitudes are 
0,15,30,45,60° from left to right; inversions are MDI OLA (top), MDI RLS 
(middle) and GONG RLS (bottom). 
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Figure 27: Rotation rates at selected latitudes and depths as a function of time, 
after subtraction of a temporal mean. The results are from GONG RLS (black) , 
MDI RLS (red), and MDI OLA (blue) inversions. 
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Figure 28: Phase (left) and amplitude (right) of 11-year sine functions fitted 
to temporal variation of the rotation rate for OLA (top) and RLS (middle) 
inversions of around 11 years of MDI observations and for RLS inversions of 
GONG data (bottom). 
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Figure 29: Local h elioseismic inferences of zonal flows close to the surface, from 
iHaber et d] (|2002h (left) and IZhao and Kosoviche^ (|2004l ) (right), reproduced 
by permission of the AAS. 
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Figure 30: Zonal flows since 1986, from Mount Wilson Doppler measurements 
(top), global helioseismic measurements from BBSO and MDI (middle) and 
MDI ring-diagram analysis (bottom). The color scale is in nHz. 
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ional flow away from the activity belts associate d with the zonal flow bands. 
Chou and Dai (|200ll) :ICh ou and Ladenkovl using data from the Taiwan 

Oscillations Network [T ON]), also found diverging m eridional flows associated 
with the activity belts. I Zhao and Kosovichev (2004) measured the zonal (Fig- 
ure [21] right) and meridional flows with the time-distance technique, and re- 
ported meridional flow converging on the activity belts above a depth of 12Mm, 
with diverging flows below 18Mm, forming circulation cells around the activ- 
ity belts. The presence of inflows i n to the activity belts was also obser ved at 
the surface bv lKomm 6^0!! (|l993bl) : iKomml ^M). iKomm et al\ (|2005l ) stud- 
ied the flows in about a year of high-resolution GONG ( "GONG-I-" ) data, and 
concluded that the overall flow pattern existed whether or not active regions 
were included in the analysis; in other words, the zonal flow bands and their 
associated converging/diverging meridional flows appear to exist independently 
of the flows in t he imm ediate vicinity of strong active regions. 

Howe et all (|2006ah compared the results from ring-diagram analysis of the 



MDI data, global analysis of MDI and GONG data, and the Mount Wilson 
Doppler observations. They found very similar results for the north-south sym- 
metrized flow pattern close to the surface in all three observations. Both the 
global and local helioseismic data indicated that the strength of the flow pattern 
did not fall off steeply below the surface. 

It should be noted that the local helioseismic observations are somewhat 
prone to systematic errors, some of which follow the changing Bq angle, or tilt 
of the solar rotati on axis relative to the observer, as shown for example by 
IZaatri et This can result, for example, in a pronounced and almost 

certainly non-solar north-south variation of the zonal flow measurements, which 
is generally corrected for by subtracting suitable averages. 

Some further features of the torsional oscillation pattern as we know it from 
a full cycle of observations from GONG and MDI (and nearly two cycles of 
surface Doppler observations) are worth noting. 

1. The exact appearance of the pattern is quite sensitive to the background 
term that is subtracted. For example, compare the /-mode results shown 
in Figure \24\ which were plotted as the difference from a smooth 3-term 
expansion of the rotation rate, with the plots in Figure 1251 which were 
plotted by subtracting the temporal mean at each location. 

2. Although the pattern repeats - of course not precisely - with each (ap- 
proximately) eleven-year activity cycle, each equatorward-migrating flow 
band exists for about eighteen years, emerging at mid-latitudes soon after 
the maximum of one cycle and flnally disappearing at the equator a cou- 
ple of years after the minimum of the following cycle; thus, the band of 
faster rotation associated with the activity of cycle 22 was still visible at 
the beginning of GONG and MDI observations in early cycle 23, and the 
band that is expected to accompany cycle 24 became visible around 2002 
(if we look at the mean-subtracted residuals), or 2005-2006 (if we use the 
smooth- function subtraction). On the other hand, each poleward-moving 
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branch seems to last only about nine years, appearing a year or so after 
solar minimum and moving to the pole before the next minimum. 

Although the equatorward-migrating bands of faster rotation are clearly 
associated with the migrating activity belts of the magnetic butterfly 
diagram, the relationship is not completely straightforward. The new 
equatorward-propagating branch is clearly visible some years before no- 
ticeable new cycle active regions begin to erupt, and the phase/latitude 
profiles of the magnetic index and th e veloc ity are very different. Also , 
as was noted bv lLabonte and Howard! 11982") and bv lHowe et all ( 2006al) . 



the strength of the torsional oscillation signal has not shown much change 
over the last few solar cycles, while the level of magnetic activity varies 
much more from one cycle to another. 

Although the equatorward branch of the zonal flow migration pattern 
shows some relationship to the pattern o f enhanced activity in the Fe xiv 



corona going back to 1973 (jAltrocla . 119971 ) . the "extended solar cycle" seen 



in these observations starts at a much higher latitude, apparently about 
70°, before migrating to the equator over about eighteen years; thus even 
the equatorward edge of these coronal activity bands seems to be at higher 
latitude than the observed new branch in the zonal flows that starts at 
about the same time. 

5. Finally, we note that because the angular velocity changes associated with 
the torsional oscillation signal are relatively small compared to the differ- 
ence in angular velocity between the surface and the bottom of the near- 
surface shear layer, while the amplitude of the signal does not decrease 
rapidly with depth, the magnitude of the shear at a given location varies 
by only a fraction of its value during the solar cycle. However, the frac- 
tional change in the shear is much greater than the fractional change in 
the rotation rate. 



9.4 Models of the torsional oscillation 



While observers, for example [Howard and Labontd (|l980 h and ' UlrichI (|200l[ ) 
have speculated that the torsional oscillation pattern might itself be part of the 
driving mechanism for the solar cycle, perhaps generating activity by shearing 
magnetic loops, modelers have generally seen it rather as a side-effect of the 

ma gnetic field s. 

Schuesslei ( 1981 ) and Yoshimural ( 198lh modeled the torsional oscillation as 



a result of the Lorentz force due to dynamo waves; according to the latter paper, 
the phenomenon would be important only close t o the surface, and wou l d have 
only equatorw ard, not pol eward, moving bands. iLabonte and Howard! (jl982[ ) 
objected to the I Yoshimural model on the grounds that it would predict a strong 
correlation between the strength of the surface magnetic field and that of the 
velocity signal, which did not seem to be the case in the observations. 
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Kueker et al. (Il996h used a different mechanism to generate the torsional 



oscillation signal in their model, considering it as the response of the Reynolds 
stress on the time-dependent dynamo magnetic field rather than a direct effect 
of the large-scale Lorentz force. This model gave a very weak poleward branch 
for the torsional oscillation signal. 

Onc e the flows had been shown observationally to penetrate well below the 
surface, iDurnev suggested that, "the pattern of torsional oscillations ap- 

pear to have the potential of critically discriminating between different dynamo 
mo dels as, e.g., the B abcock-Leighton and interface models." 

Covas et al. 1 (I2OOOI) used a model in which the observed rotation profile was 



imposed and the rotation variations arises from the action of the Lorentz force 
of the dynamo-generated magnetic field on the angular velocity. They were able 
to simulate approximately solar-like patterns of zonal flow bands and magnetic 
activity. In subsequent papers they focused on the the possibility of so-called 
"spatio-temporal fragmentation" allowing cycles of different periods in different 
regions, and in calculations wi th no density stratifi cation in the convection zone 
they found this to be feasible (jCovas et alj 2001a). The effec t was not too sen- 
sitive to uncertainties in the rotation law ( Covas et al. , 2001bl 2002 ), and some- 
what sensitive to the boundary condi tions at the o uter s urface ( Tavakol et all 
2002| ). Adding density stratification ( Covas et al . 2004 ) did not substantially 



change the results, though the amplitude of the oscillations in the deeper layers 
of the convection zone did decrease as the density gradient increased. However, 
they did find that introducing quite a small amount of a-quenching (magnetic 
fee dback on turb ulent convection) would suppress the torsional oscillation effect. 
Spruitl (|2003h modeled the torsional oscillation pattern as a "geostrophic 



flow" driven by temperature variations near the surface associated with mag- 
netic activity, and therefore having its greatest amplitude at the surface and 
falling to 1/3 of its surface value at 0.92 i?0. This model also accounts for the 
observed inflows into the activity belts. There are some problems in reconciling 
this model with the observations; it is difficult to see how the observed depth- 
dependent phase pattern could arise from a surface-originated cause, and the 
existence of the flows e ven at e pochs where there are no active regions is also 
hard to explain, though I Spruitl suggested that the flows might be produced by 
unobserved small -scale and short-lived magnetic regions. 

Rempell ( 2007 ) used a mean-field flux-transport dynamo model, with a model- 



derived differential rotation proflle and meridional flow, to investigate the ef- 
fects of various driving mechanisms for the torsional oscillation. The author 
concluded that the poleward-propagating branch of the pattern could be ex- 
plained by a periodic forcing at mid-latitudes without any underlying migra- 
tion of buried polar field. On the other hand, in this type of model the ob- 
served equatorward-propagating branch coul d not b e repro duced wit hout adding 



a therm al forcing after the manner of the ISpruitI (|2003f ) model. iHowe et al 
(20063) compared such a model with the observations, and found it not to be 
completely consistent with the observed interior behavior of the flows at lower 
latitudes. 
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10 Tachocline Variations 
10.1 The 1.3 year signal 
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Figure 31: Rotation-rate residuals at the equator at 0.72 i?0 (top) and 0.63 i?0 
(bottom), for RLS (filled) and OLA (open) inversions of MDI (red triangles) 
and GONG (black circles) data. 



Howe et al\ (|2000b( ) reported finding variations of the equatorial rotation rate 
close to the tachocline with a 1.3 year period during the early years (1995 - 1999) 
of GONG and MDI observations. The strongest signal was seen at 0.72 Rq, with 
a weaker anticorrelated signal below the tachocline at 0.63 i?©. At higher lat- 
itudes, there was also an apparent 1-year periodicity. The signal was more 
clearly seen in the GONG data, and due to the different temporal sample of the 
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MDI data it was difficult to make a quantitative comparison, but the visual ap- 
pearance of similar variations in both data sets was quite persuasive. Figure [31] 
extends the data up to the present for the equatorial locations just above and 
below the tachocline. 

Because of the role of the tachocline region in the dynamo, as well as the coin- 
cidence of the period with that seen in some heliospheric and geomagnetic obser- 
vatio ns ([Silverman and Shapiro! . 119831 : iRichardson et al. I. [l994 [Paularena et all 
19951). this claim a ttracted considerable interest, inspiring modelers such as 



Covas et al\ ([2001a[ ) to try to build models in which diffe rent periods could exist 
at t he top and bottom of th e convection zone. However. [Antia and Basu ([2000[ ) 
and [Basu and Antial ([200 1[ ). with a slightly different analysis of the same MDI 
and G ONG data, reported finding no significant variatio ns. ([Basu and Antia 
(2001) did see a signal somewhat similar to that reported bv [Howe et al. ( 2000b[ ) 
but did not consider it significant.) 

Moreover, the periodic signal disappe a rs in t he post-2001 data even in the 



original authors' analysis ([Toomre et all [2003[ : [Howe et all [20071 ). as shown 
in Figure [311 and it seems likely that the high-latitude 1-year period was an 
artifact. Intermittency in sh ort-period variations is a known phenomenon in 



the geomagnetic- index data, ([Silverman and ShapiroL [l983), and does not in 
itself imply that the phenomenon was not real. It will be interesting to see 
whether the oscillation will reappear in the new solar cycle. 



10.2 Tachocline jets 

Christensen-Dalsgaard et all (2004) searched for evidence of jet s close to the 



tachoc line, which are predicted, for example, by the model of [Dikpati et 



([2004 . Using GONG data they reported finding possible evidence of a jet at 
the tachocline, migrating equator-wards by about thirty degrees in two years 
but not at the same latitude as the surface activity belts. The significance and 
meaning of this finding remain unestablished. 



10.3 Angular Momentum Variations 

Given estimates of both density and rotation as functions of depth and latitude, 
one can calculate the solar angular momentum locally or globally. Of course, 
such calculations will reflect, and in some cases enhance, any errors in the input 
da ta, and should therefore be approached with caution. 

Komm et al ] (2003) investigated the angul a.r momentum variat ion based on 



the inversions of GONG and MDI data used bv [Howe et al. I ([2000b[ lah and found 
variations reflecting the torsional oscillation well into the convection zone and 
1.3 year variations close to the tachocline. Because the density increases steeply 
with decreasing radius, variations at greater depths will be more strongly seen in 
the angular momentum than in the rotation rate, but it should be remembered 
that no new information has been added to the data. 



Lanza[ ([2007[ ) approached the problem from the other direction, consider- 



ing the role of angular momentum transport in the modeling of the torsional 
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Figure 32: Sine- wave power in the rotation rate residuals from RLS inversions of 
GONG data, at 0.72 Rq, 0°, plotted as a function of frequency for a) 1995-2000, 
b) 1995-2003, c) 1995-2005, d) 2000-2005 
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osc illation. 

Antia et al. I (|2008l) investigated temporal variations of the solar kinetic en- 



ergy, angular momentum and higher-order gravitational multipole moments as 
derived from helioseismic inferences of the internal rotation rate; they found 
variations on the time scale of the solar cycle (but not the 1.3 year cycle), with 
some discrepancies between MDI and GONG results. They also speculate that 
the kinetic-energy changes might contribute to the observed irradiance varia- 
tions during the solar cycle; however, it is not clear that such a contribution 
is needed, as the usual view is that the solar-cycle variation in irradiance can 
be modeled simply from the effects o f sunsp ots and plage on the surface, as 
discussed, for example, bv I Jones et al\ ( 20081 ). 



64 



11 Summary and Discussion 



Since the 1970s, helioseismology has provided several insights into the interior 
solar rotation: the approximately-rigid rotation of the radiative interior; the 
differential rotation throughout the convection zone; the thin tachocline; the 
extension of the surface torsional oscillation throughout the convection zone. 
More than once, these discoveries have overtiirned theoretical expectations, in- 
spiring modelers to improve their calculations in an effort to reproduce the 
observed behavior. Because of the surprising nature of many of the findings, it 
has been important to have more than one source of observations, so that it is 
possible to distinguish between real solar features - especially the unexpected 
ones - and systematic error. 

It may be that in the future solar cycle 23, with MDI and GONG operating 
in parallel, will be seen as a golden age of helioseismology. At the time of 
writing, we eagerly anticipate the launch of the Solar Dynamics Observatory 
[SDO] with its Helioseismic and Magnetic Imager [HMI], a successor to MDI 
that will provide near-continuous helioseismic observations at higher resolutions 
than ever before and may help in unraveling the relationships between active 
region flows, magnetic fields, and geoeffective solar activity as well as providing a 
continued watch on the longer-term variations in the solar velocity fields. Sadly, 
however, current plans call for both GONG and MDI to cease to collect data 
soon after the successful launch of SDO, which would leave HMI without any 
independent cross-checks, while on the low-degree front the BiSON network has 
recently lost its funding and there are no new dedicated low-degree space-based 
instruments currently scheduled. 

There are still areas - such as the strength of the near-surface shear at high 
latitudes, the rotation of the inner core;, and any inhomogeneities and changes 
in the tachocline - that remain unclear. Furthermore, a complete numerical 
model of the solar dynamo - vital for any long-term predictive capability - is 
still lacking, and helioseismic observations still have an important part to play 
in constraining such models as they develop. 
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